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Preface 

This  study  is  the  result  of  my  attempts  to  experi¬ 
mentally  locate  instabilities  in  the  ELMAX  plasma  device. 
I  was  fortunate  to  be  given  a  sophisticated  and  well  de¬ 
signed  research  tool  already  in  operation.  With  the 
exceptions  of  a  few  leak-hunts  in  the  vacuum  system  and 
the  contruction  of  a  few  probes,  no  time  was  used  in 
preparing  equipment  for  the  experiment . 

Since  preparing  the  equipment  was  not  a  part  of  the 
study,  the  equipment'  and  data  acquisition  system  are  only 
briefly  described  in  this  paper. 

My  grateful  thanks  go  to  Dr.  Gordon  Soper,  for  his 
patient  endurance  throughout  this  study,  and  to  Major 
Justin  Curtis  and  Drs,  John  Martin  and  Merrill  Andrews 
for  their  many  suggestions  and  hours  of  assistance. 

Many  thanks  also  to  Mr.  Cliff  Van  Sickle,  Mr.  Dennis 

t 

Gross jean  and  Mr.  Bill  Cadwallender ,  without  whom  the 
sophisticated  research  would  never  have  functioned. 
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Abstract 


An  experimental  search  for  instabilities  in  an  ELMAX 

plasma  device  was  made.  Plasma  generated  oscillations  of 

potential  and  density,  with  amplitudes  several  orders  of 

magnitude  above  background  noise,  were  observed  over  a. 

wide  range  of  plasma  parameters.  Measurements  of  plasma 

parameters  were  made  with  Langmuir  probe  diagnostic 

techniques.  Frequencies  from  3  kHz  to  8  kHz  were  observed 

with  phase  velocity  in  the  direction  of  the  electron 

% 

diamagnetic  drift.  Ion  temperature  -was  not  measured  and 
attempted  measurements  of  space  potential  were  not  satis¬ 
factory.  By  comparing  the  observed  charactertistics  of 
the  oscillations  with  previously  identified  Instabilities, 
the  instabilities  responsible  for  the  plasma  generated  oscill 
ations  were  tentatively  identified  as  collisional  drift  waves 


JIL  ■ 


GEP/PH/72-22 


I.  Introduction 

The  Electrodynamic  Magnetoplasma  Experiment  (ELMAX)  is 
designed  to  achieve  improved  plasma  containment  in  a  magnetic 
mirror  device.  To  reduce  microinstabilities,  an  external 
electric  field  is  applied,  and  careful  design  of  the  plasma 
infection  devices  is  intended  to  provide  a  nearly  Maxwellian 
particle  distribution.  The  remaining  instabilities  are  ex¬ 
pected  to  be  macroscopic  in  nature  and  therefore  susceptible 
to  feedback  stabilization.  In  order  to  explore  the  feas- 
ability  of  feedback  stabilization,  it  was  first  necessary  to 
determine  what  instabilities  were  present  in  the  plasma,  and 
that  was  the  purpose  of  this  study. 

Langmuir-probe  diagnostic  techniques  were  used  to  study 
the  plasma.  The  probe  outputs  were  used  as  inputs  to  a 
spectrum  analyzer,  an  oscilloscope,  'and  a  computer  analysis 
system  developed  by  Nunn  (Ref  1).  A  description  of  the 
apparatus  is  given  in  Section  II,  and  the  experimental  pro¬ 
cedures  used,  are  described  in  Section  III. 

Since  probe  theory  in  a  magnetic  field  is  still  not  com¬ 
pletely  developed,  the  results  of  the  computer  analysis  are 
open  to  some  question.  However,  the  results  proved  to  be 
very  repeatable  and  throughout  this  study  they  are  assumed 
to  be  correct.  The  one  exception  to  this  is  the  measured 
space  potential.  As  explained  in  Section  III,  attempts  to 
use  the  gradient  of  space  potential  to  determine  the  electric 
field  gave  ambiguous  results.  Another  important  parameter 
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which  was  not  available  during  this  study  was  the  ion  tem¬ 
perature.  Although  an  ion  probe  had  been  designed  and  was 
being  built,  it  was  not  available  during  this  study.  The 
lack  of  these  two  important  parameters  made  identification 
of  the  observed  instabilities  more  difficult  and  less  certain. 
In  particular ,  it  was  found  that  precise  knowledge  of  the 
electric  field  would  permit  one  to  distinguish  the  differ- 
ences  between  a  flute  instability  p-id  a  drift  instability. 

This  was  particularly  significant  since  the  observed 
instabilities  exhibited  characteristics  common  to  both. 

Another  limitation  imposed  by  the  apparatus  was  the 

\ 

frequency  range  over  which  measurements  were  made.  Because 
of  high-frequency  attenuation  in  the  data-collecting  circuits, 
no  attempt  was  made  to  observe  frequencies  above  100  kHz, 

This  limitation  is  not  a  severe  one  since  the  most  useful 
frequencies  for  initial  experiments  with  feedback  stabili¬ 
zation  are  lower  than  100  kHz.  . 

Section  IV  includes  a  comparison  of  experimental  results 
with  theory.  If  the  drift  caused  by  the  electric  field  is 
assumed  to  be  negligible,  it  is  shown  that  all  of  the  ob¬ 
served  characteristics  of  the  instability  conform  to  those 
of  a  drift  wave,  and  that  the  stability  conditions  for  a 
drift  wave  correctly  predict  the  stability  of  the  observed 
oscillations. 
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II .  Experimental  Apparatus 

Vacuum  System 

The  plasma  was  produced  in  a  cylindrical  vacuum  chamber 
with  an  internal  diameter  of  23.5  in.  and  an  internal 
length  of  22.68  in.  The  chamber  had  three  7'* 8  in. -diameter 
ports  in  the  median  plane  spaced  90  degrees,  apart  with  one 
centered  on  the  top  of  the  chamber  and  two  on  the  sides.  One 
side  port  was  used  with  a  quartz  window  to  allow  visual 
observation  of  the  plasma.  There  were  18  radial  ports  near 
each  end  of  the  chamber,  spaced  at  20  degree  intervals, 
approximately  8.5  in.  oh  either  side  of  the  median  plane. 

The  chamber  and  other  important  components  of  the  experiment" 1 
apparatus  are  sketched  in  Pig.  1.  The  vacuum  system  included 
a  10-in.  CVC  diffusion  pump,  a  liquid-nitrogen-cooled  baffle, 
one  50  cfm  and  one  12  ofm  roughing  pump,  and  pneumatic  valves 

-7 

to  control  the  sy  :  v 'M ..  Base  pressures  of  2x10  torr  are 
obtained  using  the  baffle  and  <6x10”^  torr  without  the  baffle. 

Magnets 

As  shown  in  Fig.  1,  two  magnet  coils  were  mounted  at  each 
end  of  the  chamber.  The  two  inboard  magnets  are  4-in. -thick, 
248  turn  Magnion  coils,  They  are  mounted  next  to  the  end  of 
the  chamber  flange  with  their  centers  38.60  ±  .05  cm  from  the 
median  plane.  The  two  outboard  magnets  were  5-in. -thick,  122- 
turn  Magnion  coils  and  were  mounted  with  their  centers 
63.10  t  .05  cm  from  the  median  plane.  The  geometric  centers  of 
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the  coils  were  aligned  to  within  ±  0.05  cm  of  the  chamber 
axis.  The  power  supplies  were  able  to  provide  250  amps  of 
current  to  the  outboard  magnets  which  were  wired-  in  series, 
and  200  amps  total  current  to  the  inboard  magnets  which 
were  wired  in  parallel.  This  provided  a  magnetic  field  in 
the  chamber  of  approximately  50  gauss  at  the  chamber  axis 
near  the  end  probes  (8.5  in.  axially  from  the  median  plane) 
and  approximately  235  gauss  on  axis  at  the  median  plane. 

The  maximum  magnetic  field,  occuring  at  each  end  of  the 
chamber,  is  approximately  four  times  the  field  at  the  median 
plane,  resulting  in  a  mirror  ratio  of  /four  to  one,  At  the 
plasma  pressures  observed  throughout  this  study,  the  ratio 
of  plasma  pressure  to  magnetic  pressure  is  very  low  (Beta 
much  less  than  one)  and  the  applied  magnetic  field  is  assumed 
to  be  unchanged  by  the  presence  of  the  plasma. 

Llsitano  Coils  and  Spiral  Cathodes 

The  plasma  sources,  located  at  each  end  of  the  chamber 
between  the  inboard  and  outboard  magnets,  consisted  of  a  h— in . 
diameter  helical  Lisitano  coil  with  a  3.5  in.  i.  d.  quarts 
liner,  and  a  2  in.  diameter  heated  spiral.  The  microwave 
source  for  each  Lisitano  coil  was  a  1  kw,  2. ^5  GHs  air-cooled 
magnetron  r.f.  generator.  Both  magnetrons  were  powered  by 
the  same  ^1  kw,  500mA  d.c.  power  supply.  The  magnetrons  were 
connected  to  the  coils  through  a  wave  guide  and  stub  tuner. 
Each  Lisitano  coil  consisted  of  a  hollow  metal  cylinder  with 
a  continuous  slot  machined  through  the  cylinder  wall  in  the 
shape  of  a  helix.  The  slot  served  as  a  transmission  line  for 
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the  power  from  the  magnetrons,  and  the  net  effect  was  an 

oscillating  electric  field  applied  to  the  plasma.  The 

direction  of  the  electric  field  was  made  to  be  almost  parallel 

to  the  applied  magnetic  field  by  using  coils  with  small  helix 

angles.  The  d.c.  power  supplies  for  the  spiral  cathodes 

supplied  approximately  28  volts  at  90  amps  to  the  spirals. 

The  polarity  of  the  voltage  could  be  reversed,  and  the  spirals 

could  be  operated  with  either  end  grounded  or  with  both  ends 

\ 

floating.  Changing  the  grounding  of  the  spirals  was  found  to 
have  no  useful  effect  on  the  plasma  but  reversing  the  polarity 
of  the  voltage  resulted  in  a  reversal  of  the  applied  radial 
electric  field.  In  the  experiments  described  in  this  study, 
the  voltage  is  applied  to  the  spirals  so  that  the  radial 
field  is  directed  inward  toward  the  center  of  the  plasma.  An 
electric  field  pointing  outward  was  found  to  result  in  a 
hollow  plasma  with  maximum  density  near  the  outer  edge  of  the 
spirals  rather  than  on  the  axis  of  the  chamber. 


Langmul r  Probes 

Two  radially  traversing,  Langmuir,  single  probes  at  right 
angles  to  one  another  were  mounted  in  the  median  plane  of  the 
chamber.  Two  other  probes  were  mounted  in  two  of  the 
radial  end  ports,  one  on  top  of  the  chamber  in  line  with  t he 
top  center  probe,  and  the  second  in  a  side  port,  *10  degrees 
from  the  first.  Occasionally,  one  of  the  end  probes  was  moved 
to  the  top  end  port  at  the  other  end  of  th>  mamber ’ to  make 


measurements .  This  ailoweu  the  symmetry  of  uh  plasma  along 
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the  axis  of  the  chamber  to  be  measured,  and  also  allowed  simul¬ 
taneous  observation  of  any  observed  oscillations  at  both  ends 
and  at  the  center  of  the  chamber.  This  provided  information 
about  the  volume  of  the  plasma  affected  by  the  oscillation  and 
the  phase  relation  of  the  oscillation  at  one  end  of  the  chamber 
relative  to  the  center  or  other  end  of  the  cnamber;  The  ports 
in  which  probes  were  installed  are  marked  in  'Fig.  1  with  a  "p". 

Data  Collecting  Equipment 

To  observe  and  record  data  directly  from  the  various 
probes,  a  dual-beam  Tektronix  oscilloscope,  type  551,  was  used 
with  two  plug-in  units,  a  Nelson  Ross  spectrum  analyzer  with  a 
range  from  0  to  100  kHz,  and  a  four  channel,  chopped-beam, 
Tektronix  type  1A^  amplifier.  A  Tektronix  oscilloscope  camera 
was  used  to  record  the  data. 

Pressure  in  the  chamber  was  measured  with  an  MKS  baratron 
•pressure  meter,  type  lM.  The  flow  of  argon  gas  into  the 
chamber  was  measured  with  flow  meters  calibrated  in  standard 
'  cm  per  sec  for  argon. 

To  analyze  the  data  from  the  probes  a  computerized  system 
was  used.  The  basic  components  were  a  Hewlett  Packard  2116c 
computer,  an  HP  digital  voltage  source  (DVS),  and  an  HP  digi¬ 
tal  voltmeter  (DVM) ,  types  6130B  and  2116C.  The  DVS  and  DVM 
were  both  controlled  directly  by  the  computer.  The  DVS  was 
used  to  impress  voltage  on  a  selected  probe  and  the  amplified 
current  was  then  measured  and  fed  directly  into  the  computer. 
The  program  used  is  listed  in  Appendix  A,  with  a  brief  dis¬ 
cussion  of  the  main  features  of  the  program. 
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III.  Experimental  Procedures  and  Results 

General  Approach 

The  first  step  in  studying  plasma  instabilities  in-  the 
ELMAX  was  obviously  to  find  some  instabilities,  that  is,  to 
determine  what  values  of  the  various  control  parameters  avail¬ 
able  would  result  in  an  instability  being  observed.  An 
integral  part  of  this  step  "'as  deciding  how  to  observe  the 
instability  if  it  did  occur.  After  the  first  step  had  been 
accomplished,  the  next  step  was  to  determine  what  character¬ 
istics  of  the  instability  could  be  observed,  and  what  plasma 
parameters  could  be  measured,,  which  would  allow  the  instabil¬ 
ity  to  be  compared  to  known  types  of  instabilities  and  thus  to 
identify  it. 

Control  Parameters  -  Finding  the  Instabilities 

There  were  four  control  parameters  which  were  system¬ 
atically  varied  during  the  search  for  instabilities.  These 
were:  magnet  current,  gas  flow  rate,  current  to  the  magnetrons, 
and  current  to  the  spiral  cathodes.  In  addition  to  these 
control  parameters,  the  polarity  of  the  voltage  across  the 
spirals  could  be  reversed,  the  spiral  grounding  could  be 
changed,  and  the  magnet  current  could  be  reversed.  Variation 
of  these  control  parameters  cannot  be  directly  related  to  the 
resulting  variation  in  plasma  parameters  because  of  the  very 
complicated  interrelations  which  are  present.  For  example, 
variation  of  the  magnet  current  not  only  changed  the  magnetic 
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field  but  also  changed  the  density  and  temperature  profiles  of 

the  plasma.  In  the  same  way,  variation  of  the  current  of  the 

spirals  was  found  to  affect  not  only  the  observed  values  of 

temperature  and  density,  but  also  the  frequency  and  radial 

position  of  the  observed  oscillations.  For  this  reason  the 

control  parameters  and  plasma  parameters  were  treated  as 

distinct  from  one  another.  Control  parameters  were  varied  to 

change  the  plasma  being  studied  and  then  plasma  parameters 

measured  as  if  they  were  independent  variables.  To  detect  any 

instabilities  which  might  occur,  the  floating  potential  from  at 

least  two  probes,  usually  the  two  probes  at  the  center  of  the 

* 

chamber,  was  displayed  on  the  oscilloscope.  The  oscilloscope 
was  monitored  at  all  times  when  control  parameters  were  being 
varied.  If  an  oscillation  appeared  which  might  indicate  the 
presence  of  an  instability,  pictures  of  the  scope  trace  and 
of  the  spectrum  analyzer  display  were  made,  and  all  control 
parameters  were  recorded.  At  first,  no  useful  results  were 
obtained  using  this  technique.  Many  alternate  methods  of  * 
detecting  instabilities  were  onsidered  but  no  acceptable 
alternative  was  found.  Fc  ^unately,  during  this  time,  the 
ELMAX  system  was  being  gra  ually  improved.  Analysis  of  data 
from  previous  experiments  showed  that  higher  magnetic  fields 
than  those  being  used  in  the  present  experiment  had-  resulted 
in  a  relatively  quiet  plasma  with  what  appeared  to  be  discrete 
oscillations  sometimes  present.  This  indicated  that  increasing 
the  current  available  to  the  magnets  might  improve  results. 

The  current  available  was  increased  by  about  twenty  percent. 
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It  was  also  noted  that  the  power  supplies  to  both  the  magne¬ 
trons  and  the  spiral  had  significant  amounts  of  ripple  in  their 
output  voltage  which  might  have  been  obscuring  the  results  of 
an  instability.  Filters  .were  installed  which  significantly 
reduced  these  ripples.  With  these  improvements,  useful  oscill¬ 
ations  began  to  be  detected. 

Plasma  Parameter s-Recording  Instability  Characteristics 

Having  detected  useful  oscillations  the  next  step  was  to 
record  as  many  plasma  parameters  as  possible  while  the  instab¬ 
ility  was  present,,  and  to  observe  and  record  whatever  other 
characteristics  could  be  detected.  The  computer  analysis  of 
the  probe  curve  provided  the  means  of  measuring  the  plasma 
parameters.  The  parameters  which  were  available  from  the 
computer  analysis  included:  electron  temperature,  ion  density, 
floating  and  space  potential,  and  Debye  length.  By  varying 
the  probe  position  by  a  small  amount,  usually  one-quarter  or 
one-half  cm,  and  repeating  the  computer  analysis,  all  of  these 
parameters  were  available  as  functions  of  radial  probe  posi¬ 
tion.  This  in  turn  allows  the  calculation  of  the  gradients  of 
these  parameters.  A  typical  computer  output  sheet  is  shown  in 
Appendix  A. 

The  frequency  of  the  oscillation  was  measured  with  the 
spectrum  analyzer,  and  phase  relations  between  the  floating 
potential  oscillations  at  the  various  probe  positions  were 
recorded,  using  pictures  of  the  scope  trace.  By  biasing  the 
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probe  with  the  proper  voltage,  it  was  driven  into  the  ion- 
saturation  portion  of  the  probe  curve.  Oscillations  in  the 
current  drawn  by  the  probe  are  therefore  primarily  related  to 
ion  density  fluctuations  and  allow  a  measurement  of  the  phase 
relation  between  density  and  potential  at  that  probe.  This 
phase  relation  is  an  important  characteristic  of  the  oscilla¬ 
tion  since  different  types  of  instabilities  will  cause 
different  phase  relations  between  density  and  potential.  For 
example,  drift  wave  instabilities  usually  have  the  density 
wave  leading  the  potential  wave,  while  flute  instabilities 

normally  result  in  oscillations  with  the  opposite  density- 

\  » 

potential  phase  relation. 

The  potential  measured  in  this  experiment  is  the  floating 
potential  of  the  probe.  Since  the  density-potential  relation¬ 
ships  discussed  above  are  actually  derived  theoretically  using 
the  space  potential  of  the  plasma,  it  would  be  very  helpful  to 
determine  the  exact  phase  difference,  if  any,  between  the  space 
and  floating  potential.  This  phase  difference  is  often  * 
assumed  to  be  negligible,  and  in  the  paper  used  to  compare 
experimental  results  from  this  study  with  theoretical  predic¬ 
tions  no  distinction  is  made  between  space  and  floating 
potential.  However,  an  experiment* was  devised  which  would 
measure  the  phase  relation  between  density  .and.  floating  poten¬ 
tial  and  between  density  and  space  potential  to  insure  the 
validity  of  the  assumption  that  no  significant  phase  difference 
existed  between  the  two  potentials.  To  accomplish  these 
measurements,  a  "Boxcar"  integrator  .was  used.  Although  this 
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device  is  rather  complicated,  the  principles  of  operation  are 
straightforward.  The  boxcar  accepts  two  inputs,  the  first 
input  being  the  varying  amplitude  of  a  given  wave,  and  the 
second  being  a  trigger  pulse  synchronized  in  some  way  with 
the  input  v/ave.  The  "Boxcar”  then  samples  the  wave-  amplitude 
during  some  small  time  interval  at  a  phase  poir.t  of  the  in¬ 
coming  wave  determined  by  the  trigger  pulse.  Adjusting  the 

t 

trigger  pulse  to  coincide  with  consecutive  phase  points  along 
the  incoming  wave  allows  the  characteristics  of  the  wave  to  be 
determined  as  a  function  of  phase  position.  In  this  experiment, 
the  input  wave  was  the  output  of  the  Langmuir  probe  as  the 
biasing  potential  v/as  varied  to  generate  the  probe  curve.  The 
output  of  the  "Boxcar"  was  then  used  as  the  input  to  the 
computer  which  analyzed  the  probe  curve.  In  this  way.  the 
various  plasma  parameters,  specifically  the  temperature,  density, 
and  potential  were  calculated  as  a  function  of  phase  position  of 
the  oscillation.  By  plotting  these  variables  versus  the  phase 
position,  the  relation  between  these  variables  could  be 
determined.  The  results  of  this  experiment  were  inconclusive 
for  two  reasons.  First,  several  hours  were  required  to  complete 
one  run  through  one  wavelength  of  the  oscillation.  When  the 
beginning  of  the  run  v/as  repeated  Immediately  after  completing 
a  run,  to  check  the  repeatability  of  the  results,  different 
results  were  obtained.  Second,  a  careful  analysis  of  the  output 

of  the  "Boxcar"  revealed  that  the  output  was  not  maintaining 
a  constant  relation  to  the  input  amplitude  as  it  should  have. 
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The  results  of  this  particular  experiment  we re  therefore  not 
used. 


Experimental  Results 

Continuous,  almost  sinusoidal  oscillations  were  observed 

over  a  wide  range  of  control  parameters.  These  oscillations 

had  an  apparent  clockwise  phase  velocity  as  looking  along  the 

magnetic  field  and  in  the  direction  of  the  magnetic  field. 

To  facilitate  discussion,  the  direction  of  the  magnetic  field, 

which  was  parallel  to  the  longitudinal  axis  of  the  vacuum 

chamber,  is  designated  as  the  positive  z-direction.  Then  the 

\  • 

r  and  theta  directions  are  designated  as  the  other  two  direc¬ 
tions,  the  apparent  phase  velocity  of  the  oscillations  as  in 
the  positive  theta  direction.  Wave-length  along  the  z-axis 
could  not  be  measured  exactly,  since  probes  could  be  inserted 
at  'only  three  points  along  the  z-axis.  If  the  three  probes 
along  the  z-axis  were  all  adjusted  to  the  same  radial  position, 
the  center  probe  showed  the  largest  amplitude;  however,  if  the 
three  probes  were  adjusted  to  the  radial  distance  which  gave 
the-  largest  amplitude  for  that  probe,  the  three  amplitudes  were 
approximately  the  same.  For  the  oscillations  which  are 
analyzed  in  detail  in  this  study,  the  center  probe  had  the 
largest  oscillation  at  a  radial  position  of  3  cm.  The  two  end 
probes  showed  the  maximum  oscillations  when  positioned  at  2.5  crn. 
When  the  probes  were  adjusted  to  these  positions,  no  phase 
difference  between  the  three  probes  could  be  detected.  These 


observations  can  be  exolained  if  a  longitudinal  scanuinc  wave 


is  assumed  to  exist  in  the  plasma.  The  boundary  conditions  for 
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this  standing  wave  would  be  the  physical'  characteristics  of  the 
chamber.  Since  the  first  conducting  surface  which  was  in  con- 
•  tact~wi~th  the  plasma' at  each-  end  of  the  'chamber  was  fhe '‘spiral"”’ 
cathode,  the  wavelength  of  the  standing  wave  must  be  some  inte- 
gral  multiple  of  the  distance  between  the  two  spirals .  A 
wavelength  equal  to  this  length  would  result  in  a  minimum 
amplitude  at  the  center  of  the  chamber.  Since  this  was  not 
observed,  the  wavelength  of  the  standing  wave  must  equal  at 
least  twice  the  distance  between  the  spirals. 

To  determine  the  wavelength  of  the  oscillations  in  the 
theta  direction,  the  output  of  two  probes  in  the  same  r-theta 
plane  were  compared  on  the  oscilloscope,  using  a  copped  sweep 
on  the  oscilloscope  so  that  the  phase  relation  between  the  two 
signals  could  be  seen.  When  using  the  center  probes,  which  were 
positioned  90  degrees  apart,  a  phase  relation  of  either  90  or 
180  degrees  was  observed.  The  change  from  one  phase  relation  to 
the  other  was  always  accompanied  by  an  apparent  doubling  of 
frequency  and  was  caused  by  a  change  in  the  applied  magnetic 
field.  This  was  interpreted  as  a  mode  change.  The  term  "mode" 
as  used  here  refers  to  the  ratio  between  the  circumference  of 
the  plasma  column  at  the  radial  position  at  which  the  oscilla¬ 
tion  is  detected  and  the  wavelength  of  the  oscillation.  Since 
the  oscillation  is  observed  to  be  continuous,  the  wavelength 
must  be  equal  to  the  circumference  divided  by  some  integral 
number.  If  this  were  not  true,  interference  would  occur  at 
some  point  around  the  plasma  column  and  a  continuous  oscillation 
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would  not  bo  observed.  If  Xj>  is  the  wavelength  of  the  oscilla¬ 
tion  in  the  theta  direction,  and  Kj^  is  the  wave  number,  then 

K:  =  27r(A|)"1  =  2ti  ( 27TP/M)**1  (1) 

K|  =  2TTM/2Trr  =  M/r  (2) 

where  M  is  the  mode  number.  When  an  apparent  phase  shift  of 
90  degrees  is  observed  betv/een  two  probes  positioned  90  degrees 
apart,  this  is  interpreted  to  mean  that  the  same  maximum  poten¬ 
tial  arrives  at  the  second  probe  one-fourth  of  a  cycle  after  it 
reached  the  first  probe.  This  is  called  an  m=l  mode.  If  the 
phase  shift  appears  to  be  180  degrees,  this  means  that  one  probe 
is  sensing  a  minimum  potential  while  the  other  is  sensing  a  maxi 
mum.  This  is  interpreted  to  mean  that  the  wavelength  of  the 
oscillation  is  equal  to  one-half  of  the  circumference,  that  is, 
m=2. 

The  sketches  in  Fig.  2  illustrate  the  results  of  a  sequence 
of  spectrum  analyses  of  the  oscillation.  The  input  to  the 
spectrum  analyser  was  the  floating  potential  from  a  probe  at 
the  center  of  the  chamber,  radially  positioned  to  the  point  of 
maximum  amplitude  of  the  oscillatipn.  The  magnetic  field  was 
increased  after  each  analysis  was  made  so  the  sequence  illus¬ 
trates  the  behavior  of  the  oscillation  as  a  function  of  magne¬ 
tic  field.  Fig.  2(a)  shows  a  mode  1  oscillaf^n  with  a 
frequency  of  S  kHz.  With  the  magnetic  fie"1  i ;  r eased,  Fig  2(b) 
shows  a  decrease  in  the  frequency  of  the  -  >;*  _  oscillation  to 
6  kHz  and  at  1?  kHz  car.  be  seen  the  beginning  of  a  mode  2 
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Pis.  2  Spectrum  Analyses  of  Oscillations 
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oscillation.  Fig.  2(c)  shows  the  continuing  decrease  in  fre¬ 
quency  of  both  modes  as  the  magnetic  field  is  increased.  The 
increase  in  the  amplitude  of  the  mode  2  oscillation  and  the 
decrease  in  the  amplitude  of  the  mode  1  can  clearly  be  seen 
in  Fig.  2(d)  and  2(e).  In  Fig.  2(e)  the  mode  1  oscillation 
has  almost  completely  disappeared,  while  the  mode  2  oscillation 
has  reached  its  maximum  amplitude  and  continued  to  decrease  in 
frequency . 

To  determine  the  phase  relation  between  the  density  and 

potential  oscillations,  a  probe  was  biased  with  a  negative 

voltage  so  that  the  current  drawn  by  the  probe  was  a  measure  of 

\ 

the  ion  density  in  the  plasma  surrounding  the  probe.  The 
current  was  then  amplified  and  displayed  on  the  oscilloscope  to 
determine  the  phase  relation  between  density  and  potential.  The 
density  was  found  to  lead  the  potential  by  approximately  90 
degrees,  both  in  mode  1  and  mode  2. 

The  oscillations  were  observed  over  a  wide  range  of  control 
parameters.  The  limits  were  usually  imposed  by  equipment  limit¬ 
ations  rather  than  the  loss  of  the  oscillations.  The  magnetic 
field  ranged  from  200  gauss  to  250  gauss  in  the  center  of  the 
chamber.  No  oscillations  were  observed  below  200  gauss  and  250 
gauss  required  the  maximum  magnet  purrent  available.  Oscilla¬ 
tions  were  found  only  in  the  upper  twenty  percent  of  the  range 
available  in  spiral  current.'  This  range  corresponds  to  the 
largest  density  gradients  which  were  found  in  plasma  density. 

Pressure  could  be  varied  from  about  0.^  microns  of  H  to  about 

g 

2.5  microns  without  losing  the  oscillations,  but  were  primarily 
studied  at  a  pressure  between  1.5  to  2.0  microns  because 
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the  conditions  were  the  easiest  to  maintain  in  this  pressure 
•region.  Control  of  the  power  to  the  magnetrons,  and  adjustment 
of  the  stub  tuners  was  not  used  primarily  to  vary  the  plasma 
parameters,  but  rather  were  both  adjusted  so  as  to  insure  equal 
production  of  plasma  at  both  ends  of  the  chamber  with'  a  minimum 
amount  of  noise  in  the  magnetron  circuits.  This  was  done  to 
minimize  the  presence  of  source-generatpd  oscillations  in  the 
plasma. 

Reversing  the  current  to  the  magnets  resulted  in  no  apparent 
change  in  the  oscillations,  other  than  a  reversal  of  the  theta' 
phase  velocity  relative  to  the  chamber,  as  expected. 

Reversing  the  polarity  of  the  spiral  voltage  resulted  in  a 
hollow  plasma  with  maximum  density  occuring  at  about  3  cm  rather 
than  at  the  center  of  the  plasma.  This  plasma  configuration  was 
not  studied  because  of  the  large  amount  of  noise  present  in  the 
probe  outputs  and  because  of  the  undesirable  density  profile. 

With  this  knowlege  of  the  general  characteristics  of  the 
oscillations,  a  general,  qualitative  comparison  of  the  instability 
responsible  for  the  oscillations  could  now  be  made.  For  a 
quantitative  comparison,  however,  some  information  on  the  vari¬ 
ation  of  the  plasma  parameters  wab  required.  To  obtain  this 
information,  radial  scans  of  the  plasma  electron  temperature, 
the  ion  density,  and  the  space  potential  of  the  plasma  were  made. 

This  was  done  by  making  a  computer  run  with  the  probe  in  a 
given  radial  position  and  obtaining  a  computer  output  which 
gave  the  desired  parameters  at  that  probe  position.  _  The  probe 
was  then  moved  to  a  position  a  small  distance  from  the  first 
position  and  another  computer  run  .was  made.  The  plots  shown  in 
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Fig.  3  through  Fig.  8  display  the  ion  density  and  electron  tem¬ 
perature  as  a  function  of  probe  position.  Fig.  3  shows  the 
parameters  as  measured  at  the  top  probe  with  a  mode  1  oscilla¬ 
tion  present  in  the  plasma.  Fig.  *1  gives  the  parameters  with 
the  magnetic  field  reduced  until  the  oscillations  disappear. 

Fig.  5  gives  the  same  information  as  Fig.  1 1  except  that  the 
oscillations  were  removed  by  reducing  spiral  current  rather 
than  magnetic  field.  Fig.  6  through  Fig.  8  give  the  same  inform¬ 
ation,  except  that  the  end  probe  was  used  rather  than  the  center 
probe. 

The  significance  of  the  data  given  in  Fig.  3,  ll ,  and  lj  can 
.best  be  seen  by  comparing  the  data  observed  with  the  oscilla¬ 
tions  present  to  each  of  the  two  cases  with  the  oscillations 
not  present.  The  temperature  variation  is  very  slight  in  either 
case,  however  the  density  with  oscillations  present  (Fig.  3)  is 
three  times  as  great  as  with  no  oscillations  due  to  reduced  mag¬ 
netic  field  (Fig.  4)  but  only  half  again  as  large  as  the  density 
when  the  oscillations  are  removed  by  reducing  the  current  to  the 
spiral  cathodes.  As  is  discussed  in  comparing  the  observed 
oscillations  with  known  types  of  oscillations,  the  density 
gradient  emerges  as  the  important  factor  rather  than  the  density 
itself.  The  same  relations  are  seen  in  Fig.  6,  7,  and  8  for  the 
data  taken  fr.in  the  end  probes. 

Although  space  potential  is  calculated  in  the  program  used, 
it  is  not  plotted  because  it  was  found  not  to  be  useful  in  deter¬ 
mining  the  nature  of  the  instability  causing  rl  served 
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oscillations.  Near  the  end  of  this  study,  some  improvements 
in  the  computer  program  was  made  which  caused  different  values 
to  be  calculated  for  space  potential.  This  occurred  too  late 
for  new  runs  to  be  made,  however,  so  no  valid  information  about 
the  space  potential  is  available  for  this  study. 

Another  plasma  parameter  which  v/ould  aid  in  a  quantitative 
comparison  of  the  observed  instabilities  with  known  theory  is 
the  ion  temperature.  No  method  wa s  available  during  this  study 
to  measure  ion  temperature. 

Th  frequency  observed  ranged  from  3  kHz  to  8  kHz  and 
exhibited  a  strong  dependence  on  magnetic  field.  As  the  magne¬ 
tic  field  is  increased,  the  observed  frequency  decreases  until 
a  mode  change  occurs.  At  this  point,  the  frequency  doubles  and 
then  begins  to  decrease  again  with  increasing  magnetic  field. 
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Pig.  7  End  Probe  -  Reduced  Magnetic  Field 
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Fig.  8  End  Probe  -  Reduced  Spiral  Current 
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IV ;  ‘  Theory 

The  very  large  number  of  known  instabilities  makes  a 
comparison  of  a  given  experimentally  observed  instability 
with  all  known  instabilities  almost  impossible.  The  litera¬ 
ture  survey  by  Auer,  et  al*.  proved  invaluable  in  making  the 
comparison.  This  survey  lists  136  known  instabilities  and 
cites  references  for  each.  Many  of  the  instabilities  listed 
could  be  eliminated  immediately  because  of  some  obvious  dis¬ 
agreement  with  the  observed  characteristics.  For  example, 

all  of  the  ion-cyclotron  instabilities  were  eliminated  just 

\ 

on  the  basis  of  the  frequency  observed  being  much  too  low  for 
this  type  of  instability.  (Ref  2) 

Many  other  instabilities  would  exhibit  experimental 
characteristics  which  are  in  agreement  with  those  observed, 
but  can  be  eliminated  on  the  basis  of  some  one  calculated 
plasma  parameter,  such  as  one  of  the  collision  frequencies. 

An  example  of  this  type  is  the  Simon-Hoh  instability  (Ref  • 
3:383)  which  was  eliminated  because  the  ion-neutral  collision 
frequency  is  too  low  in  the  experimental  plasma  (Ref  ^:14). 

The  Kelvin-Helmholtz  instability  was  considered  to  be 
almost  identical  to  the  observed  Instability  until  the  parallel 
wavelength  was  measured.  Since  the  Kelvin-Helmholts  wave  is 
characterized  by  short  parallel  wavelengths,  the  observation 
of  long  parallel  wavelengths  allowed  elimination  of  this 
instability  (Ref  5:818). 
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Temperature-gradient-driven  drift  instabilities  also  ex¬ 
hibit  characteristics  very  similar  to  the  observed  waves, 
however,  they  are  stabilized  by  the  presence  of  a  density 
gradient  and  one  stability  criterion  is  a  simple  relation 
between  the  density  and  temperature  gradients-  which  the  exper¬ 
imental  plasma  does  not  satisfy  (Ref  6:1107). 

There  were  two  types  of  instabilities  which  were  consid¬ 
ered  which  could  not  be  eliminated.  These  were  the  collisional 
drift  instability  and  the  flute  instability.  The  term 
"collisional"  as  used  here  refers  only  to  the  mechanism  which 
restricts  the  motion  of  electrons  parallel  to  the  applied 
magnetic  field.  As  a  convenient  measure,  the  drift  instability 
is  considered  collisional  if  the  mean  free  path  of  the  electrons 
between  collisions  with  ions  is  less  than  the  parallel  wave¬ 
length  of  the  oscillation.  Flute  Instabilities  are  almost 
always  considered  probable  in  a  mirror-type  geometry  unless 
some  provision  such  as  Ioffe  bars  are  employed  to  suppress 
them.  Drift  instabilities  have  been  found  in  almost  all  .kinds 
of  confinement  schemes,  and  are  in  fact  so  common  that  they  are 
often  called  "Universal  Instabilities."  The  two  instabilities 
are  very  similar  to  each  other,  and  in  fact  the  flute  is  often 
treated  as  a  special  case  of  the  drift  instabilities. 

(Ref  7:154-156) 

Quantitative  Comparison  V.’ith  Flute  Instability 

The  dispersion  relation  for  the  flute  instability  as  given 
by  Schmidt  (Ref  7:283-294)  is  used  in  this  analysis.  This 
dispersion  relation  leads  to  the  stability  condition 
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K2/ft2  (g  +  u2^1/^)2  >  i\g  n 1/n()  (3) 

where  K  is  the  perpendicular-  wave  number,  ft  is  the  ion 
cyclotron  frequency,  g  is  the  effective  gravity,  uf„  is  the 
thermal  speed  of  the  ions  and  nVn Q  is  the  normalized  den¬ 
sity  gradient.  If  the  plasma  is  assumed  to  be  nearly  thermal, 
and  no  radial  electric  field  is  considered,  then  g  =  (^)  U^/R. 
If  the  inequality  in  Egn.  3  is  not  satisfied,  unstable  oscilla¬ 
tions  result .  Using  data  from  the  top  center  probe  where  flute 
instabilities  are  most  likely  to  occur  and  assuming  ion  temper¬ 
ature  to  be  about  300  degrees  Kelvin,  the  values  for  the 
experimental  plasma  are 

K  =  1  x  io3  m"1 

5.5x  10^  rad  sec"1 

p  li  p  „p 

U*  =  6.2  x  ioH  m  sec  * 

li  _p 

g  =  1.9  *  10  m  sec 
n1/nQ  =  ^9  m"1 

and  the  inequality  is  not  satisfied.  However,  solving  the 
dispersion  relation  for  the  real  part  of  the  frequency  gives 
a  value  of  less  than  one  kHz  which  does  not  agree  with  the 
observed  frequency.  In  addition,  since  the  drift  velocity 
Imparted  to  the  electrons  and  ions  by  a  flute  instability  is 
proportional  to  the  mass  of  the  particles,  the  drift  of  the 
electrons  is  negligible  compared  to  the  ion  drift  velocity  and 
therefore,  the  oscillation  travels  in  the  direction  of  the  ion 
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drift  velocity.  This  results  in  the  potential  wave  leading 
the  density  wave  which  is  the  reverse  of  the  observed  relation. 
Both  of  these  difficulties  can  be  overcome  if  a  radial  electric 
field  is  assumed  such  that  the  plasma  column  rotation  would 
give  an  apparent  observed  frequency  of  4  kHz.  With  such  an 
electric  field,  the  doppier  shift  due  to  the  plasma  E*B  drift 
would  dominate  the  oscillation's  phase  velocity  and  the  density- 
potential  phase  relation  would,  be  reversed.  In  such  a  case 
however,  the  value  of  g  would  be  changed  by  the  addition  of 
a  term  V2/R  where  V  is  the  tangential  velocity  of  the  plasma 
and  R  is  the  distance  from  the  center  of  rotation.  To  give 
the  proper  frequency,  V2/R  =  9.6xl0^m  sec"2  and  g  becomes 
9.78xio^rn  sec“2.  When  this  new  value  for  g  is  inserted  into 
Eqn.  3,  the  inequality  is  satisfied  and  no  instability  will 
occur.  Justification  for  the  added  term  for  the  effective 
gravity  is  given  by  Bhatia  (Ref  9:1652).  The  resulting 
stability  is  predicted  by  Vlasov  (Ref  10:486). 

* 

Quantitative  Comparison  With  Drift  Instability 

During  the  literature  search  for  this  study,  a  paper 
written  by  Hendel  and  Politzer  (Ref  11)  was  found  which 
appeared  to  give  almost  identical  results  to  the  ones  ob¬ 
tained  during  this  study.  Subsequently,  a  later  paper  by 
Hendel,  Chu  and  Poliozer  was  found  which  derived  a  disper¬ 
sion  relation  explaining  the  results  of  the  earlier  paper. 

The  observed  characteristics  of  wave  number  (both  parallel 
and  perpendicular),  frequency  range,  density-potential  phase 
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relation,  and  the  dependence  of  frequency  and  mode  number  on 
magnetic  field  we  re  all  in  agreement  with  .the  observations 
made  by  Hendel,  et .  al .  (Ref  12)'. 

The  dispersion  relation  used  by  Hendel,  et  al.  is  de¬ 
rived  using  linear  theory  and  a  slab  model  for  the  plasma. 

In  the  regime  of  parameters  found  in  the  experimental  plasma 

in  this  study,  the  resulting  stability  condition  reduces  to 

1 

K11  (nl/n0)“2  >  (me/mi)2  (4) 

where  is  the  parallel  wave  number,  nVn^  is  the 

normalized  density  gradient,  m  is  the  electron  mass  and 

a 

mi  is  the  mass  of  the  ions.  Using  data  from  the  end  probes, 
this  stability  condition  not  only  predicts  .instability  when 
the  oscillations  are  observed,  but  also  predicts  stability  when 
the  oscillations  are  eliminated  by  reducing  either  spiral 
current  or  magnetic  field.  The  correct  predictions  are  con¬ 
sidered  the  best  indication  that  the  observed  oscillations  are 
indeed  caused  by  a  collisional  drift  instability. 


For  a  numerical  analysis,  Eq(^)  can  be  rearranged  to 


read 


K11  (me/tni)  2  >  (nl/n0)2  (5! 

1 
2 

For  argon  atoms,  (me/m.  >  equals  3.69*10  .  If  X ^  is 

p 

equal  to  twice  the  distance  between  the  spiral  cathodes, 

_  jj  _  p 

equals  6.85*10  cm  ,  and  the  stability  condition  becomes 


0.185  cm"2  >  (n3/n0)2  cm"2 
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or  the  equivalent  statement,  the  wave  grows  if 

(n^/n^J  >  0JJ30  cm“^  (7) 

Since  the  maximum  amplitude  oscillations  were  observed 
at  a  radial  position  of  3.0  cm  for  the  center  probes  and 
2.5  cm  for  the  end  probes,  ( n1/nQ |  is  calculated  at  these 
positions  from  the  data  given  in  Pigs.  3  through  8.  For  the 
center  probes  Jn1^)  equals: 

0.^9  -  with  oscillations 

0.12  -  with  reduced  magnetic  field 
a$d  no  oscillations 

0.28  -  with  reduced  spiral  current 
and  no  oscillations 

For  the  end  probes  | n1/nQ |  equals: 

0.80  -  with  oscillations 

0.23  -  with  reduced  magnetic  field 
and  no  oscillations 

0,20  -  with  reduced  spiral  current 
and  no  oscillations 

The  theory  agrees  with  observation  in  all  cases  and  oscilla 
tions  are  observed  only  when  the  inequality  in  Eon.  (7)  is 
satisfied. 
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V.  Conclusions  and  Recommendations 

The  purpose  of  this  study  was  to  find  experimentally  and 
to  idencify,  using  existing  theory,  the  dominant  instabil¬ 
ities  in  the  ELMAX  plasma  device.  The  good  agreement 
between  the  theory  developed  by  Hendel,  et.  al..  was  very 
satisfying.  The  recommendations  which  are  suggested  by  the 
results  obtained  are  simple  and  obvious.  The  need  to  be  able 
to  measure  ion  temperature  and  parallel  wave-length  were  in 
fact  already  known  and  under  consideration  before  the  study 
began.  The  study  did  point  out  the  need  for  a  better  way  to 

V 

measure  the  space  potential  so  that  the  electric  Held  in  the 
plasma  could  be  calculated. 

The  follow-on  experiment  is  already  under  way  and  some 
initial  success  with  suppression  of  the  observed  oscillations 
lend  some  support  to  the  conclusions  reached.  In  particular, 
suppression  of  the  oscillations  using  suppressor  probes 
located  near  the  end  of  the  chamber,  where  the  observed  density 
gradient  is  the  largest  and  where,  therefore,  the  oscillations 
are-  predicted  to  be  most  likely  to  grow,  lends  support  to  the 
conclusions  reached. 

The  real  conclusion  of  this  study  must  wait  until  the 
experiments  on  feedback  stabilisation  have  been  completed. 
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Appendix  A 

Computer  Program  Used  In  Analysis 

The  computer  program  used  in  the  calculation  of  electron 
temperature  and  ion  density  is  listed  In  FORTRAN  in  this 
appendix, • immediately  following  fig.  9.  Pig.  9  shows  a 
typical  computer  output  from  the  program.  The  program  is 
designed  to  analyze  the  characteristics  of  a  probe  curve, 
which  is  the  plot  of  probe  current  versus  probe  voltage.  The 
program  is  written  to  use  a  variable  number  of  points  in  '  ,e 
analysis  but  normally  250  points  were  used.  In  Pig.  9  the 
first  numerical  entry  is  simply  entering  the  number  of  poir'"' 
to  be  used.  To  generate  the  probe  curve,  a  minimum  and  maximum 
voltage  are  selected  by  trial  and  error  to  include  both  ion 
saturation  current  and  electron  saturation  current.  In  the  run 
which  resulted  in  Pig.  9>  these  values  are  -23  volts  and  +1 
volt.  The  next  entries  tell  the  computer  which  probe  is  being 
used  and  the  radial  position  of  the  probe.  The  last  entry 
refers  to  the  gain  on  an  amplifier  used  in  the  probe  current 
circuit . 

Since  the  probe  current  is  the  sum  of  the  electron  and  ion 
current,  a  straight  line  approximation  of  the  ion  current  is 
made,  and  then  subtracted  from  the  total  current  to  give 
electron  current.  The  approximation  used  is  the  next  item  given 
in  Pig.  9.  Next  follows  a  listing  of  selected  currents  and 
probe  voltages  with  the  "INDEX"  showing  at  which  of  the  250 
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points  they  occured.  This  read-out  helps  in  the  initial  sel¬ 
ection  of  the  maximum  and  minimum  voltages,  used  and  is  also 
helpful  in  analysing  the  operation  of  the  computer. 

To  calculate  the  electron  temperature,  a  non-dimensional 
potential  is  defined  and  the  electron  current  equation  is 
written  as  an  exponential  function  of  this  potential.  Since 
the  defined  potential  is  a  function  of  electron  temperature, 
a  plot  of  the  natural  log  of  electron  current  versus  potential 
yields  a  straight  line,  the  slope  of  which  determines  the 
electron  temperature. 

Since  the  current  at  ion  saturation  is  a  function  of  the 

\  « 

ion  density,  a  value  for  density  can  be  calculated  once  the 
ion  current  is  known.  In  this  program  the  density  is  calcula¬ 
ted  several  ways  and  four  different  densities  are  shown  on  the 
computer  output.  In  this  study  the  one  selected  is  calculated 
using  theory  developed  by  Laframboise  assuming  ion  temperature 
to  be  zero  (Ref  13),  The  selection  of  this  particular  density 
is  not  too  important  since  in  this  study  the  absolute  valuo  of 
the  density  is  not  important,  but  rather  the  density  gradient. 
All  of  the  densities  given  by  the  program  yield  approximately 
the  same  gradients .  The  electron  temperature  is  shown  in  Fig. 

9  as  "TE  at  ( I'FL) "  and  is  given  in.  degrees  .Kelvin.  The  density 
appears  as  the  last  entry  under  "DENSITY." 
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ENTER  RUN  NUMBER 

1 

enter;  number  of  points 

VMIN, VMAX 

PROBE  NUMBER, POSIT I ON  (CM) 
CONVERTER  GAIN 

250 

-23,1 

2,-1 

.1 


ION 

CURRENT 

APPROXIMATION 

( 

51.9  INCREMENTS) 

Cl 

=  .1593 

E+0 1  *  V  +  - . 

4557E+02  E 

MAX  =  .86G0E+01 

ZERO  ION  CURRENT  AT  28. 

60 

VOLTS 

INDEX 

'  V 

CP 

CES 

MIN 

1 

-• 1533E+02 

. 6445E+02 

MI 

106 

- . 5205E+0 1 

- 

•5675E+02 

MM1 

113 

-•4530S+0 1 

•» 

•  5 1 20E+02 

•3101E+01 

MM3 

120 

- • 3855E+0 1 

- 

•498OE+02 

•3196S+01 

IFL 

160 

. 000OE+0O 

. 1200E+01 

.6176E*02 

ISP 

203 

•4 145E+0 1 

• 2260S+04 

. 2446E+G4 

MM4 

236 

• 7325E+01 

• 5225E+04 

•4669E+04 

MM2 

243 

•8000E+01 

.5321E+04 

•4858S+04 

MAX 

250 

•8675E+01 

• 5594E+04 

TE  AT  V(IFL) 

=  12515.  DEG 

K 

EMAX  = 

•4339E-0 1 

SHEATH  CALCULATIONS 

14.07  VOLTS  BELOW  V(ISP>  (146.0  INCREMENTS) 

BETA2  =  . 1978E+01 

SHEATH  THICKNESS  =  .7032E-O4 

RADIUS  =  . U45E-C3  METERS 

LENGTH  n  .9786S-03  METERS 

RP/RS  =  *  3659E+0S 

ION  NUMBER  DENSITY  = 

•  124 19E+13  AT  V(IFL) 

•50325E+ 12  AT  V(IFL)-10  (USING  SHEATH) 

DEB  TI/TE  RP/LD  VPF  DENSITY 
•869E-03  1  • 5 12E+0 1  4.93  .789E+12 

. 773S-03  0  • 574E+0 1  5.24  .999E+12 

PAUSE 
PAUSE 


PiC.  9 


Typical  Compute:'  Output 
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oeoi 

0038 

C003 

oso4 

0005 

0006 

0007 

0308 

0009 

0010 

0011 

0012 

0013 

COM 

0015 

0016 

0317 

0018 

0019 

0020 

C021 

0022 

0023 

0024 

0025 

0026 

0027 

0023 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0033 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

OMSO 


ftn.b  Appendix  A 

PROGRAM  DAT A9 

COMMON  VC25O>,C23CO5O),C2C250>.AC3>,BCS>.CC6>»ALC2SC25O> 
*TSC250>.SDC2SC2S J).DC2>.THC2>.TXLC2) 

■  CALL  CLPS 
MUL=3 

E=1.6021E-19 
XM=  1.33052-23 
TEN6= 1 .2-G6 
AC1>=0. 15113 
AC  2) =0.14661 
AC3>=0. 13303 
AC4>=C. 1 1041 
AC5>=G. 073733 
AC6>=0. 038009 
AC7)«-D«01 1765 
AC 8 >*-0.070583 
DC  1 >  =  0.0 
DC2)=3. 57142-3 
BC3>*7. 14292-3 
BC4>*1. 07142-2 
BC5)*1. 42862-2 
BC6>»1. 78572-2 
BC 7 >=2.14292-2 
BC8>*fr.025 
CC 1 >=-9.04932-3 
CC2>=-8. 56502-3 
CC3>=-7. 11052-3 
CC4> =-4.68652-3 
CC5)=-1. 29232-3 
CC6>=3. 07052-3 
CC7> =3.40342-3 
CCB>=1. 47062-2 
XK*40. 

XK!^  1.67252-27 
XME*9. 10912-31 
PI =3. 14159 
EPS=8, 55532-12 
VKITEC  2.  1 000> 

1000  FORMAT  C7/"EMT2H)  DA7E"/7X. "DIAMETER  AND  LENGTH  CIN  IN 
♦CUES)  OK  PROBE  1"/7X. "DIAMETER  AND  LENGTH  OF  PROBE  2'V 
♦7X."VIS .VKP.VSP") 

READC  1.  103ONUM1  A.NUKIB 
1030  K0KMATC2I3> 

DO  180  1*1.2 
READ  Cl.*>  TRCO.TXLCI) 

TRC I >=TRC I >*2.542-2/2. 

180  TXLCI >=rXLC I >*2.542-2 
READC l.*>V IS. VFP.VSP 
CALL  L2ADUC8.20> 

10  PAUSE  1 

WRITEC2, 1010> 

1010  FORMAT  CS/"2>J7SR  RUM  NUMBER'* > 

READC i.*>MUM2 
40  ELITE  C2,1020> 


1020  FORMAT  C/"2>!TSP.;  MUMPER  OK  POINTS,,/7X.',VMIM.VMA:c"/7X. 


♦••PROBE  NUMBER. POSITION  C  CM  >"/7.X,  "CONVERTER  GAIN") 


READ  Cl.*>  ME 


READ  Cl.*>  VC1>,VCMU 
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0060 

HEAD  <1,0  IP, CM 

0061 

READ  <1,0  GAIN 

0062 

R=TJ;<  IP) 

0063 

xl=txl<ip) 

0064 

DV*<V<MD-V<  1  >  >/FL0AT<NWd  > 

0065 

DO  20  I =2, NO 

0066 

20 

V<  I  )=V< 1-1 )+DV 

0067 

vs=v<N'.:>-v<i> 

0068 

CALL  SLV.  <14,  ID) 

0069 

IF  < ID)  31,32 

0070 

31 

1/1:1  TS  <2,7100) 

0071 

7100 

FORMAT  < //"ANALYTIC  GSNERAT ION"/"ENTER  XI, X2  FOR  RANGE 

0072 

♦OF  X") 

■0073 

RE  AD<  1,0X1, M2 

0074 

DX=<X2-X1  >  /FLOAT  <NW) 

0075 

X»X1 

0076 

DO  710  Id, Mi: 

0077 

CPEL=5000«  *<  1 •+TANHCX)  ) 

0078 

CPI0N=5.*X-1C. 

0079 

CPHIE=20**<  1 .  +T AMH ( 3 • *  <  X+2 •  >  >  ) 

0030 

CE<  I  )=CPSL+CPIOM+CF:!IS 

0031 

710 

X»X*DX  . 

0032 

DO  25  1=1, MM 

0083 

25 

CE< 1 )=CS( 1 )/GAIN 

0034 

GO  TO  5 

0035 

32 

CONTINUE  ^ 

0086 

DO  33  Id,NK 

0087 

33 

CS<I)=0.  •  / 

0088 

0089 

GAIMdO.*GAlN 

MC0DE=4 

WHITE  <  10)  NCODS  SLAVS' 

DO  30  I d,Nw 

0090 

0091 

0092 

CALL  DVS  <V<I>> 

0093 

DO  35  Jd,lO 

0094 

HEAD  <10,1040)  IS,TCE 

0095 

1040 

FORMAT  <  IX,  11,310*0) 

0096 

ISdS+1 

00.97 

GO  TO  <35,2,3, 3), IS 

0098 

2 

TCS«-TCS 

0099 

GO  TO  35 

0100 

3 

WRITE  <2,1041)  I 

0101 

1041 

F05-.MAT  </"DVK  OVERLOAD  AT  DATA  POINT",  14) 

0102 

CALL  DVS  <V<1)> 

0103 

GO  TO  40 

0104 

35 

CE< 1 >=CS< I >+TCS 

0105 

IF  <10.-AnS<CE<I>)>  1,4 

0106 

1 

NC0DE=3 

0107 

WHITS  <10>  NCODS 

0103 

4 

CE< I )«CE< I )/GAIN 

0109 

CALL  HOPS  <  IEF.it) 

0110 

IF  <ISRH-l)  30,8 

0111 

8 

WRITS  <2,1045)  I 

0112 

1045 

FORMAT  </"DVS  ERROR  AT  DATA  POINT", 14) 

0113 

GO  TO  40 

0114 

30 

CONTINUE 

0115 

CALL  DVS  <V<1>) 

0116 

S 

CALL  SWR<0, ID) 

0117 

IF  < ID)  7,9 

0118 

7 

CALL  PPLOT  <V,CE,NW> 

0119 

CALL  3NR<0,ID> 
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0120 

IF  <ID>  40,9 

0121 

9 

IFL»1 

0122 

DO  60  I»2,NV; 

0103 

IF  <AD3<CE<I)>-Ar>S<CE(IFL>>)  11,60  . 

012/1 

11 

IFL=I 

0125 

60 

CONTINUE 

0126 

cpf=cecifl> 

0127 

VPFl=V< IFL) 

0128 

CALL  DUS  (VPF-1  > 

0129 

DO  70  i*i,*v: 

0130 

70 

vmsv<i>-v?Fi 

0131 

DVIS=VIS/DV 

0132 

IDV«=DVIS 

0133 

IF  (IDV-3)  12,13 

013/1 

12 

IDU=3 

0135 

13 

IF  ( IFL-< 1BV+1 ) >  14,15 

0136 

14 

WRITE  {2,  1C60> 

0137 

1060 

FORMAT  •< "V(MIM)  TO  U!S  INCLUDES  V(IFL>") 

0138 

GO  TO  10 

0139 

IS 

CALL  t.ILSQ  CV<1 >,CS< 1 >, IDV,D,ER1,SK2> 

0140 

SL«D(1> 

0141 

CP0=D<2) 

0142 

VPCs-CPO/SL 

0143 

WHITE  <2,  1070)  DVIS,SL,CP0,EP.t,EH2,VPC 

0144 

1070 

FOi.MAT  CURRENT  APPROXIMATION  <",F5.1,"  INCREXEN 

0145 

*tsi»/2:<,"ci  io. 4,"  *•  j  +",eii.4,/«:<,»emax  =",sio.4. 

0146 

*3:<,"Kr.::s  «",s  10.4m,  "zero  iom  cup.pjcmt  at»,f6.2,»  volts" 

0147 

IKA5C»IFL-1 

0148 

DO  80  I=IDV, I  MAX 

0144 

J=IFL-I 

0150 

CIC«Q»99*<CP0+SL+V< J)  ) 

0151 

IF  <ClC-CE< J) )  80,21 

0152 

80 

CONTINUE 

0153 

WHITE  <2, 10301 

0154 

1080 

FORMAT  < /"STRAIGHT  LINE  APEOX  MOT  VALID") 

0155 

GO  TO  10 

0156 

21 

MI«J 

0157 

DO  90  !*MI,NV 

0158 

IF  <V<I)-VPC)  17,15 

0159 

17 

TS< I )=C3< I )»(CPO+SL*VC I >  > 

0160 

18 

IF  (TS<I>-0.01>  19,90  - - gp% 

0161 

19 

TS<I)r.O,Ot  15# 

CONTINUE 

INONUL-1  faeP'0 

MM  1  »M  1 4 1 NC  \ 

0162 

0163 

90 

0164 

0165 

KM2=NV:-INC 

0166 

DO  110  i«=m:u,mm2 

0167 

S1«=A(  1  )*TS<  I  ) 

0168 

DO  100  J»1,INC 

0169 

K«J+l 

0170 

IP=I+J 

0171 

IM«I-J 

0172 

SUM*TS(IM)+TS(IP> 

0173 

100 

S1«=SUA<K)*SUM 

0174 

IF  <S1-C.01>  101,102 

0175 

101 

Sl=0.Ol 

0176 

102 

CEi>( I )=8 1 

0177 

110 

ALC.“S<  I  )  =ALOGCS  1 ) 

0178 

MM3«:*MW1NC 

0179 

MM4t=M:*2-INC 
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0180  DO  130  I»MM3,XM4 

0181  S2«0. 

0162  S3=C( 1 )*CES< 1 ) 

0183  DO  120  <J=  1,  IMC 

0184  K=J+1 

0185  ip=i+.j 

0186  IM=l-J 

0137  S2«S2+n<K)*(ALCE5( IP)»ALCSS( IK)> 

0133  120  S3=S3+C(K)*<CSS(  IM)+CSii(  IP) ) 

0189  S2=S2/D V 

0190  IF  <-S2)  121# 122 

0191  122  S2=1C.0.*2/XK 

0192  121  T3{  I  )«2/<S2*.TiO 

0193  S3=S3/<DV*DV) 

0194  130  SBC2S<I)=S3 

0195  MM4A«MK4-5 

0196  DO  140  t  =  IFL,K>'4A 

0197  IF  (SDCSS(  I  >  )  131,140 

0193  131  IF  (SDCES< 1+5 )  >  132,140 

0199  132  IS?**  I 

0200  IF  <SDCES( I« 1 )+SDCES( I  >  )  133,150 

0201  133  IS?=*ISP-l 

0202  GO  70  150 

0203  140  COMTIWJE 

0204  ISP=IFL 

0205  WHITE  <2, 1200) 

0206  1200  FORMAT  <"DATA  DOES  DOT  EXTEflD  70  PLASMA  POTENTIAL") 

0207  150  CORfIMUE 

0208  WHITE  <2, 1090) 

0209  1093  FORMAT  </6X,"!MPSX",7X,"V",  18X,"CPM#  IOX,"CES"> 

0210  WHITE  <2, 1100)  V<1),C2<1) 

0211  1100  FORMAT  <"MIM",6H,"1",  1X,2S13.4) 

0212  WHITS  <2,1110  1 1 1  ,V<Ml  ),CS<M1 ) 

0213  1110  FORMAT  <"MI",I3.- 1X,SE13. 4) 

0214  WHITE  (2,1120)  MM  1 , V < MM 1 ) , CE < MM  1 ) , CSS  <MM 1 ) 

0215  1120  FORMAT  ("MM l", l 7, ID, 3E13.4 > 

0216  WRITS  <2,1130  !:M3,V<MM3),CE(MM3),CSS<KM3> 

0217  1130  FORMAT  <"K>13",  17, 1X,3S13.4) 

0218  WHITS  <2,1140  IF1.,V<  IFL),CE(  IFL),CES(  IFL) 

0219  1140  FORMAT  <"IFL",I7,  1X,3S13.4> 

0220  WHITE  (2,1150)  ISP,V< ISP), CS< ISP),CES< ISP) 

0221  1150  FORMAT  (•'ISP*',  17,  1M,3E13.4 > 

0222  WRITE  <2,1160)  MM4 , V  <  MM4 ) , C S  <  MM4  > , C E S  <  MM4 ) 

0223  1160  FORMAT  <"MM4",I7, 1X,3S13.4) 

0224  WRITE  <2,1170)  MM2 , V < MM2 ) , CS < MM2 ) , CSS <KM2 ) 

0225  1170  FORMAT  <"KM2",I7, IX,3S13.4> 

0226  WHITE  <2,1130  NW,V<ML >,CS<MW> 

0227  1180  FORMAT  ("MAX", 1 7, 1X,2E13.4> 

0223  DVFP=VFP/DV 

0229  IDVbDVFP 

0230  IF  (IDV-3)  151,152 

0231  151  IDV«3 

0232  152  M5T=IFL+IDV 

0233  CALL  LILSO.  <V<MST),ALCES<MST), IDV,D,EH1,SF.2) 

0234  IF  < -D( 1 ) >  153,154 

0235  154  DU^IOO.+S/XX 

0236  153  TSS«=X/<X;<+D<  1)> 

0237  WHITE  <2,1190)  7K5,.'H1,ER2 

0233  1190  FORMAT  </"TS  AT  V<IFL)  =",F7.0,"  DEG  K",5X,"SMAX  «",E10 

0239  *.4,3X,"ERMS  *=",E10.4) 
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0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

155 

0251 

156 

0252 

0253 

171 

0254 

172 

02S5 

157 

0256 

0257 

160 

0258 

159 

0259 

0260 

0261 

0268 

0263 

0264 

162 

0265 

0266 

161 

0267 

164 

0268 

0269 

163 

0270 

165 

0271 

0278 

0273 

0274 

2001 

0275 

027S 

0277 

0278 

0279 

167 

0280 

0281 

168 

0888 

1210 

0283 

0284 

158 

0885 

0836 

0287 

0888 

0289 

0290 

501 

0291 

0292 

0293 

0294 

0295 

0296 

299 

0297 

* 

0298 

VDL*=V(ISP)*D(  i> 

Cl =*(CPO+SL*V{ JFL)  > 

DEUS2=fl. 

DEN2A=0. 

DE\’S4  =  0. 

DENSl=Cl*S0R7(XM*XMP/<XK*TES)>*l.S-l2/(1.13l‘S*PI<tR*XL) 

DVSC=10.0/DV 

!DV=DVSC 

ICB=IFL-IDV 

IF  < - ICS)  156,155 

icr>=i 


CP=CS( JCB) 

IF  (IC3-MI)  172,171 
C?=C?04-SL*V(1C3) 

IF  (CP)  157,15=5,157 
V?1C2=A5S(V(ISP>-V(ICD>) 
IF  (ISP-IFL)  159, 16ft, 159 
vpicb=v?icr+5«tt:-:s*>:x/s 


DVSC=VPIC-VDV 

DELA=SGRT(E*(2.*V?1CS)**3/(XM*X>?P)) 

DELDs4.*PI*EP5*XL/(9.*iV!‘C?vl.S-6) 

DSLCsSORT (  E*V ? !  C«V  <  XK* TES  )  ) 

BE7A2=AnS(DELA*DZL3*(  1.+8.66/DELC)) 

IF  (0.5-3STA2)  161,162 

DELS=R+(9.55472-2*3S7A2+t.00l*S0RT(nETA2)-5.7212E-5) 


GO  TO  165 

•IF  (5.0-3ZTA2)  163,164 

D2LS=R*(7. 8  1S-2*2STA2+ 1 , 03  15*S0RT (BETA2)- 1 . 383E-2) 


GO  TO  165 

DF.US=;:*<6*8S4S-2*13S7A2*l.0791*SGK7(BSTA2)-7.413Sn-2> 

RS=R+DSLS 

xs=xl+dels 


RD=R/ES 

WHITE  (2,2001  >  VPIC3,DVSC,BSTA2,DELS,RS',XS,nD 


FORMAT  (/••SHEATH 
*V(1SP)  (",F5.l," 
•SHEATH  THICK-MESS 
*"/5X, "LENGTH  =", 


CALCULAT 1 OX5"/5X , F5  *  2,  "  VOLTS  BELOW 
INCR2M2XTS)"/5X,  "BETAS  =",E  10.4/SX, " 
=",210. 4/5X, "RADIUS  =",S10.4,"  METERS 
El 0.4,"  MSTSKS'VSX, "IIP/RS  *",S10.4)< 


IF  <  14. -32TA8)  163,167 

DSNS4  =  -  ( CP*S  QET( XM*XNP/  <  XK*TES  >  >  *  1 .  E-  1 2  )  /(  1 .  1 3*E* 
*PI*i!S*X5) 

WHITE  <2, 1210)  DEM SI, DENS 4 

FORMAT  (/"ION  Mu MEEK  DENSITY  ="/5X,E11.5,"  AT  V(IFL>"/ 
*5X,E11.5,"  AT  V< IFL)“ 10  (USING  SHEATH)") 

CONTINUE 
BX=Xi(*l.E7 
RpnR+100. 

EV=XK*TES/S 
CGE=4. 00222- 10 
WHITE  (2,501) 

FORMAT  ( /3X, "DEB", 4X, "7 I /TS», 4X, "HP/LD", 5X, "VPF", 5X, 
♦"DENSITY") 

RT*S<MiT (  XME/ ( MM*XN? )  ) 

DEL'  1  =SORT  {  5I(*TE3/(  4«  *PI  *CGE*CGE*DENS  1 )  ) 

JP=0 

VLT=10.0 
JPASS=(! 

DZD=DZ31 
HPLD=EP/DSB 
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0300  300  GO  TO  <301, 303, 310), J? 

0301  301  ITITS»1 

0302  C !  I  =2. E“4*HPLD*F.PLD-R.  03 1 752*hPLD*2.8366 

0303  GO  TO  305 

0300  303  ITITSsO 

0305  Cll  - 1 . 5S-/»*.".Pr.D*:;FLD-0.  02343*fiPLD*2.257fl 

0306  305  VPM=«ALOG<KY*Cl  1 ) 

0307  vpf=vp:;*sv 

0308  IVS=IFL*IF  t:«VPr/PV) 

0309  IVKaVLT*EW/CV 

0310  IVI=1V3-IVN 

0311  VIVI«V<IVI> 

0312  IF  (-IVI)  307,306 

0313  306  VIVI=V<t)+FLOAT<IVI-l)*DV 

0314  307;  CPP=CFO+SL*VIVI 

0315  DENS3*-CPP*SS:rt'<XK*XMP/<a.*Pl*XK*TSS))*i.E-12/<CX:*R*XL*E> 

0316  DEB3=SGr;T<BK*TSS/<4.*?I*CG:!*CGK*DENS3>> 

0317  DEBD«An$<DSn-DSB3> 

0318  IF  <0200-0.01*0231  >  309,314 

0319  314  023*0233 

0320  kpld*::p/deb 

0321  JPASSaJPASfi+1 

0322  IF  <C-.J?A5S>  309,300 

0323  309  VISITS  <2, 502)  6233,ITl?E,Ii?LD,VPF,DENS3 

0324  502  FORMAT  <23. 3, IS, 212. 3, F8. 2, 21 1.3) 

0325  IF  <9-J?ASS)  311,299 

0326  311  VISITS  <2,500)  0230 

0327  500  FOl'.KiiT  <"**0200  *",F6.4> 

0328  GO  TO  299 

0329  310  CONTINUE 

0330  PAUSE  4 

0331  CALL  SVtS  <1,ID) 

0332  IF  < 10)  22,23 

0333  22  N»ISP-MM3*5 

0334  CALL  PPLOT  <V<KM3),TS<MM3),M) 

0335  23  CALL  SUI  <2, ID) 

0336  IF  < ID)  24,26 

0337  24  MBITS  <4,1050)  NUMlA,i\'UMlD,NUM2,CM,TES,VFPl,V<IVS>,DSNS3 

0333  *,D2NS4 

0339  1050  FORMAT  <313,6210.4) 

0340  26  CONTINUE 

0341  GO  TO  10 

0342  END 

0343  ENDS 


•♦END-OF-TAPE 

* 
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